Coherent inelastic backscattering of intense laser light by cold atoms 
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We present a nonperturbative treatment of coherent backscattering of intense laser light from 
cold atoms, and predict a nonvanishing backscattering signal even at very large intensities, due to 
the constructive (self-)interference of inelastically scattered photons. 

PACS numbers: 42.50.Ct, 42.25.Dd, 32.80-t, 42. 25. Hz 



When a plane wave of arbitrary nature is incident upon 
a disordered medium of scatterers, the backscattered in- 
tensity is an interference pattern of all coherent partial 
amplitudes containing detailed information on the sam- 
ple configuration. Under an ensemble average, interfer- 
ence between uncorrelated amplitudes is washed out, ex- 
cept for a small angular range around exact backscat- 
tering, where the average intensity may exhibit a nar- 
row peak. This peak results from constructive inter- 
ference between multiple scattering probability ampli- 
tudes counterpropagating along direct and reversed paths 
P, 0- This phenomenon is called Coherent Backscat- 
tering (CBS) and was for the first time demonstrated 
with samples of polystyrene particles 0. The CBS en- 
hancement factor a, the ratio of the total intensity at ex- 
act backscattering to the background intensity, measures 
the coherence of counterpropagating amplitudes respon- 
sible for localization effects. Recently, CBS of light has 
been imported to the quantum realm with clouds of cold 
atoms 0,01 El An important leitmotiv of these studies is 
the robustness of the underlying interference effect with 
respect to fundamental quantum mechanical dephasing 
mechanisms, such as spin-flip (of the incoming radia- 
tion, which carries a polarization degree of freedom) 
or inelastic scattering. This has important repercussions 
for the transition from weak to strong (in Anderson's 
sense) localization of light in disordered atomic samples 
, and also for potential technological applications such 
as random lasers . While in the regime of weak, per- 
turbative atom-field coupling, the partial destruction of 
CBS due to spin-flip like processes - induced by the 
multiple degeneracy of the atomic transition driven by 
the incident radiation - has been demonstrated exper- 
imentally and analysed theoretically in quite some de- 
tail, experiments and theory only now start to probe the 
strong coupling limit, where inelastic photon- atom scat- 
tering processes prevail. First experimental results on 
Sr (driving the Sq — > 1 Pi transition with its nondegen- 
erate ground state, hence in the absence of spin-flip) 6] 
indeed demonstrate the reduction of the CBS enhance- 
ment factor with increasing intensity of the injected field, 
for values s — f2 2 /2(A 2 + 7 2 ) < 1 of the atomic satura- 



tion parameter (where fl is the Rabi frequency induced 
by the driving, 7 half the spontaneous decay rate of the 
excited atomic level, and A the detuning of the injected 
laser frequency from the exact atomic resonance) . A first 
scattering theoretical treatment ^| identified the origin 
of such suppression in the availability of which-path infor- 
mation through inelastically scattered photons: reversed 
paths can be distinguished by the detection of photons 
of different frequency. However, this treatment, still per- 
turbative in the field intensity, cannot address the limit 
of large saturation parameters s > 1 - with an emerging 
Mollow triplet ^lj m the single atom resonance fluores- 
cence - and, in particular, makes no prediction on the 
crossover from dominantly elastic to essentially inelastic 
CBS, nor on the CBS enhancement factor in the deep in- 
elastic limit. In the present Letter, we enter this regime, 
starting from a general master equation which allows for 
a nonperturbative treatment of the atom-field coupling. 
As we will see, even inelastically scattered photons give 
rise to a nonvanishing CBS signal. 

We start out from the elementary toy model of CBS - 
a laser field scattering off two atoms with labels 1 and 2, 
placed at a fixed 01 distance ri 2 3> A = 27r/fc£, with 
the wave vector of the incident field. It is known from 
the perturbative treatment of CBS that double scattering 
(on two atoms) provides the leading contribution to the 
CBS signal, since this is the lowest order process which 
gives rise to time reversed scattering amplitudes which 
can interfere constructively. We expect that this scenario 
also allows for a qualitative assessment of the nonlinear 
atomic response in the regime of high laser intensities, 
whilst propagation effects in the bulk of the scattering 
medium are certainly beyond reach of this model. Dis- 
order will be mimicked by a suitable average over the 
atomic positions. Our model also neglects the accelera- 
tion of atoms out of resonance, which certainly becomes 
important at very high intensities, but can be experi- 
mentally compensated for by shortening the CBS probe 
duration, as realized in ||. We focus exclusively on the 
photon coupling to the internal atomic degrees of free- 
dom. 

With these premises, the average backscattered in- 
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FIG. 1: (Color online) Elementary configuration for coherent 
backscattering (CBS) of intense light (thick arrows) by two 
isotropic dipolar transitions in the helicity preserving polar- 
ization channel (dashed arrows). The sublevels |1) and |3) of 
both atoms have magnetic quantum number m = 0; sublevels 
|2) and |4) correspond to m = — 1 and m = 1, respectively. 



tensity can be derived from the correlation functions 
of the atomic dipoles which emit the detected signal 
[Til Hij | . We specialize to the scenario of the Sr exper- 
iments, with a nondegenerate atomic dipole transition 
J g = — » J e = 1, driven by laser photons with right cir- 
cular polarization on the sublevels |1) — > |4) - see Fig.^ 
The scattered light is detected in the helicity preserv- 
ing channel (i.e., of photons which originate from the 
\2) —> |1) transition), where single scattering is absent. 



Thus we obtain, up to an irrelevant prefactor, the expec- 
tation value for the stationary intensity scattered into the 
direction k close to the backward direction — k^, 

(!)„ = (<4) ss + (a| 2 ) ss + 2Re(( ( r 2 1 l( T2 2 ) ss e lk - r -) , (1) 

where <rjy = \k) a (l\ a , for atom a. The steady-state 
values for correlation functions of the form (ag) ss or 
(ofjO^i )sb can be found from the master equation 0] 



(2) 



where the Liouvillians C a and govern the evolution 
of an arbitrary atomic operator Q for independent and 
dipole-dipole interacting (through the exchange of one 
or several photons) atoms, respectively. Q stands for an 
operator from the complete set of operators acting on a 
tensor product of Hilbert spaces of individual atoms. For 
our choice of the atomic structure shown in Fig. ^ 
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The explicit form of the interaction-picture Liouvillians 
C a and La/5 derived in the standard dipole, rotating- 
wave, and Born-Markov approximations can be shown to 
read: 



C a Q = -iApt • D a , Q] - i[n o (Dt • e L ) + fi* (D Q ■ e* L ), Q]+ 1 fa a ■ [Q, DJ + [D* , Q] ■ D Q ) , (4) 
C a0 Q = Dt • T(g, n) ■ [Q, T> ] + [D+ , Q] ■ T*(g, ft) D Q , (5) 

I 



where A = ujl — ujq is the detuning, fi a = Sle lkL ' r ° is 
the atomic (coordinate-dependent) Rabi frequency, and 
£ l fixes the polarization of the laser field. 

D a = -e_i<T™ 2 + £ oo"i3 - (6) 

is the lowering dipole operator of atom a, with e±i, £o 
the unit vectors of the spherical basis. The radiative 
dipole-dipole interaction due to exchange of photons be- 
tween the atoms is described by the tensor T (g, ft) = 
jg A , where A = 1 — fift is the projector on the plane 
defined by the vector n connecting atom 1 and 2, and 

pikori2 

<? = *3rr > (7) 

2fc ri2 

where fcrj = Wo/c, is the small coupling constant \g\ -C 1 
in the far-field limit fcrjri2 3> 1, where we neglect near- 
field interaction terms of order l/(fco r i2) 2 and l/(fco.ri2) 3 



(which, at higher atomic densities, could also be retained 
in our formalism). 

Transforming the operator equation @ to a system of 
255 linear coupled differential equations for the atomic 
correlation functions, we can solve for the physical quan- 
tities which enter the expression for the detected in- 
tensity. In doing so, we furthermore take advantage of the 
far field limit fco r i2 3> 1, and expand the correlation func- 
tions up to second order, ((. . . , in the dipole-dipole 
coupling constant J7|). The double scattering contribu- 
tion to the CBS signal, detected in the helicity preserving 
channel, is then precisely given by terms proportional to 
g 2 , since it stems from the exchange of two photons be- 
tween the atoms, along a 'direct' and its 'reversed' path. 
Finally, the CBS signal is obtained after an elementary 
configuration average (. . .) CO nf. defined through the fol- 
lowing twofold procedure: (i) isotropic averaging of the 
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FIG. 2: (Color online) Total intensities of the ladder, the 
crossed terms, and their sum 7s°* ' 2 ' , in the helicity preserving 
channel, as functions of the saturation parameter s at reso- 
nance, A = 0. 



relative orientation f i 2 over the unit sphere; (ii) uniform 
averaging of the distance ri 2 over an interval of the or- 
der of A, around a mean value given by the mean free 
path. After this simple procedure all terms relevant for 
the calculation of the CBS enhancement factor survive, 
whereas all the irrelevant terms vanish. 

We thus arrive at our final expression for the total 
second-order intensity 



(8) 



a sum of the total ladder (or background), L tot , and total 
crossed (or interference) term C tot (8), with 8 the obser- 
vation angle of the scattered intensity with respect to the 
backward direction. In terms of the second order atomic 
correlation functions, L tot and C tot (9) are given by 



L tot = ((- 2 1 2 )I 2 s 1 + («W 
C to \9) = 2Re/(a 2 y i2 >I 2 s ] e lk - r " 



conf . 



(9) 
(10) 



Therefrom we deduce the main quantifier of CBS, the 
enhancement factor 



C tot (0) 



(11) 



Figure 12 shows our results for the total CBS intensity 
as well as its components L tot and C tot (0), as a function 
of the saturation parameter, at exact resonance (A = 0). 
The behavior of 7*°* ' 2 ' shows that the double scattering 
intensity behaves markedly different from that of an iso- 
lated atom. Whilst the scattering intensity from an iso- 
lated atom tx s/(l + s) is known to saturate for large 
s 1 1 5| - the double scattering intensity exhibits a maxi- 
mum at s ~ 0.7, followed by gradual decrease ~ for 
large s: At high laser intensities, more and more photons 
are scattered inelastically, and are therefore less likely to 
undergo resonant interaction with the second atom. 



FIG. 3: (Color online) Enhancement factor a in the helic- 
ity preserving channel, versus saturation parameter s. Solid 
curve: on resonance (A = 0), dashed curve: off resonance 
(A = 7). Straight lines represent the perturbative predic- 
tion 2 - (1 + S)s/4 of [[J. Inset: The finite enhancement 
Qfoo ~ 1.09 signals residual photon (self-)interference, even in 
the deep inelastic regime. 



The enhancement factor a(s) follows directly from 
the above quantities. As shown in Fig. [3J it decays 
monotonously from its weak field limit a(0) = 2. In qual- 
itative agreement with the experiment this decay is 
faster for finite detuning A = 7. For small values of s, a 
is well approximated by the linear decay 2 — (1 + 5)s/4, 
with 5 = (A/7) 2 , derived within the scattering picture 
fiof . For large values of s (inset), however, a satu- 
rates at a value ~ 1.09 strictly larger than unity, 
whilst one would expect vanishing contrast £Le., a = 1) 
for scattering from two independent atoms [17], Hence, 
the (self-)interference of inelastically scattered photons 
unambiguously contributes to the crossed term C tot (0). 
Note that this observation bears some similarity to CBS 
with photons from degenerate Raman transitions, which 
were shown to yield an important contribution to the 
CBS contrast, even in the limit of infinite ground state 
degeneracy |lj| , as well as to the residual CBS enhance- 
ment in optically active media at high magnetic fields 
0. 

In contrast, elastically scattered photons remain per- 
fectly coherent, and contribute to the CBS intensity with 
a constrast two, for any s. To see this, we just need to 
extract the purely elastic component of the signal from 
the total yield in cq. QJ. Since the detected intensity 
{I)ss is nothing but the autocorrelation function of the 
source field amplitudes radiated by the atomic dipolcs, 
its elastic part (J)^ is generated by the classical dipolcs 
induced by the injected radiation - this is by their av- 
erage, nonfluctuating parts (cr^)ss .15]. Hence, (I)f s 
is given by the product of the expectation values of the 
atomic dipoles: 



(I)f s = |<a 2 1 1 ) ss | 2 + |(a 2 2 1 } ss | 2 + 2Re(( ( T 2 1 1 } ss (a 1 2 2 ) ss e lk - r -). 

(12) 

A power series expansion of the right hand side of i|12D 
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to second order in the coupling g leaves only symmetri- 
cally factorized combinations of the form (c^i)' 1 '^^)^. 
Asymmetric combinations, like (c^i } ' 2 ' i^u ) do not 
contribute to the signal since the |1) <-» \2) transitions 
are not laser-driven (see Fig. QJ, hence (crf 2 )'°' vanishes. 
Evaluation of the correlation functions by symbolic cal- 
culus, together with the configuration average described 
above, finally provides an analytic expression for the elas- 
tic ladder and crossed terms: 

L* = C-(0) = 24^1^^. (13) 

Expression shows that the elastic ladder and 

crossed terms are equal for any s, as to be expected 
from reciprocity arguments |l3j . These elastic compo- 
nents decay like s~ 3 at large saturation, much faster than 
the total intensities that decrease like s _1 (cf. Fig. 
This proves that the residual CBS enhancement is 
entirely due to the (self-)interference of inelastically scat- 
tered photons. Furthermore, expression l|13fl shows that 
the elastic part of the double scattering intensity exhibits 
a maximum at s = 1/3, slightly below the departure of 
a(s) from the perturbative prediction of 0| m Fig- 01 
Consistently, an expansion of (|13fl to second order in s 
reproduces the expression L cl = C el (0) ~ s — 4s 2 de- 
rived in ^(J- Note that the crossover to the nonlin- 
ear regime for double scattering occurs at a value of 
s three times smaller than for an isolated atom, where 
7 cl [°] ocs/(l+s) 2 exhibits a maximum at s = 1. This 
has a transparent interpretation, by virtue of factorizing 
eq. into (i) the elastic intensity / cl [°] scattered by 
the first strongly driven atom, (ii) the total scattering 
cross section a tot tx 1/(1 + S)(l + s) of the second atom, 
and (iii) the relative weight pl[o]/jtot[o] = ^y^tot = 
( 7 2 + A 2 )/( 7 2 + r> 2 /2 + A 2 ) = 1/(1 + s) of elastic 
processes therein. Obviously, higher order scattering pro- 
cesses than considered in our present contribution must 
unavoidably push the crossover value of s to even smaller 
values. 

In conclusion, we have presented the first study of 
coherent backscattering of intense laser light from sat- 
urated dipole transitions. The CBS enhancement de- 
creases monotonously as a function of s, but, remarkably, 
coherence is - partially - preserved in the deep inelas- 
tic limit of the two-atom response to intense laser radia- 
tion, since also inelastically scattered photons can inter- 
fere with themselves, along time-reversed paths. Conse- 
quently, CBS should also have an imprint on the spec- 
trum of the scattered radiation, as well as on its photo- 
count statistics, which are both directly accessible in the 
framework of our present approach, as well as in labo- 
ratory experiments. Furthermore, let us note that our 
present results are also relevant in the somewhat differ- 
ent context of Young's double slit experiments with two 



atoms 17]. In contrast to the forward Young-type inter- 
ference that necessarily decoheres for s — > oo, since the 
photon visits two different, uncoupled atoms, a backscat- 
tering experiment, with appropriate polarization sensi- 
tive excitation and detection, must lead to a finite inter- 
ference contrast. 
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